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A. 


FUNDAMENTAL  STUDIES  OF  ROCKET  COMBUSTION  CHEMISTRY 


A.l  RESEARCH  OBJECTIVES 

The  combustion  of  solid  rocket  propellants  and  other  energetic  materials  is  a  complex  multidimensional 
and  multiphase  process  involving  a  wide  variety  of  chemical  species.  The  very  high  pressure  and 
temperature  conditions  of  practical  rocket  combustion  chambers  are  at  present  inaccessible  by  most 
conventional  diagnostic  techniques.  The  study  of  these  coupled  phenomena  in  situ,  therefore,  has  not  been 
possible  in  sufficient  detail  to  develop  a  complete  understanding  of  the  chemistry  and  physics  of  the 
combustion  process.  The  objective  of  most  recent  studies  of  the  combustion  of  these  materials  has  been 
to  study  separate  aspects  of  the  overall  process  in  an  effort  to  provide  a  comprehensive  understanding  of 
the  combustion  mechanism.  This  study  is  one  component  of  that  coordinated  investigation  and  has  as  its 
focus  the  gas  phase  reactions  associated  with  the  combustion  of  these  solid  fuels. 

The  decomposition  of  many  of  these  solid  energetic  materials  during  combustion  leads  to  the  formation 
of  gaseous  hydrocarbon  fuel  species  and  oxides  of  nitrogen  which  serve  as  oxidizers  (ref.  1  and  2).  The 
reactions  of  these  decomposition  products  above  the  propellant  surface  lead  to  a  gaseous  flame  which  can 
provide  heat  which  is  transferred  back  to  the  propellant  surface  and  can  thereby  influence  the  burning  rate. 
The  purpose  of  this  paper  is  to  summarize  the  current  status  of  studies  we  have  undertaken  of  model  gas 
phase  flames  associated  with  the  combustion  of  nitramine  bases  solid  rocket  propellants.  These  studies 
consist  of  measurements  of  the  structure  of  stable  and  unstable  species  concentration  profiles  and 
temperature  in  laminar,  premixed,  flat  flames  of  fue!/NOx  mixtures  at  low  pressure.  The  experimental 
measurements  arc  then  compared  to  calculations  of  the  concentration  profiles  using  a  one  dimensional 
flame  code  which  models  the  transport  processes  and  chemistry  of  the  flame.  The  transport  processes 
include  species  diffusion  and  thermal  conduction  through  the  flame  and  the  chemistry  is  modeled  bv  a 
detailed  chemical  kinetic  reaction  mechanism. 

The  basic  mechanism  used  for  the  modeling  work  is  a  subset  of  the  331-reaction  mechanism  o!  Volpom 
and  Branch  (ref.  3).  Reactions  removed  were  those  involving  species  with  low  concentrations,  including 
condensed  phase  molecules,  and  molecules  containing  more  than  three  carbon  atoms.  The  resulting 
mechanism  contains  272  reactions.  This  mechanism  was  the  basis  for  all  of  the  modeling  described  in 
this  paper.  The  flame  code  used  to  solve  the  one  dimensional  flame  equations  was  by  Kce  et  al.  (ref.  4) 


A.2  RESEARCH  APPROACH  AND  RECENT  RESULTS 
A.2.I  FLAMES  SUPPORTED  BY  N,0  AS  OXIDIZER 

Several  recent  studies  have  focused  on  measuring  the  temperatures  and  species  concentrations  in  flames 
supported  by  N:0.  Habeebullah  et  al.  (ref.  5),  Zabarnick  (ref.6),  and  Vandoorcn  el  al.  (ref.  7)  studied  CH . 
reacting  with  N;0.  Habeebullah  measured  stable  species  concentrations  using  probe  sampling,  radical 
species  concentrations  using  laser-induced  fluorescence  (LIF),  and  flame  temperatures  using  thermocouples 
and  LIF.  Zabarnick  used  LIF  to  measure  stable  and  radical  species  concentrations  and  flame  temperatures. 
Dindi  et  al,  (ref.  8)  studied  CO  reacting  at  N;0,  and  used  probe  sampling  to  measure  stable  species 
concentrations,  LIF  to  measure  unstable  species  concentrations,  and  thermocouples  to  measure  flame 
temperatures. 

We  have  modeled  the  structure  of  all  of  these  flames  using  our  272  step  reaction  mechanism.  It  was 
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found  that  in  order  to  have  best  agreement  with  the  full  range  of  data  represented  by  these  results,  the 
rales  of  several  of  the  reactions  needed  to  be  adjusted  as  indicated  in  Table  1.  The  most  significant  ot 
these  changes  was  for  the  N;0  decomposition  which  was  adjusted  to  within  the  upper  range  of  the 
literature  data  for  this  reaction.  This  change  was  needed  to  model  the  methane  flame  data  correctly.  The 
other  reactions  listed  in  Table  I  were  adjusted  to  provide  good  agreement  for  the  CO/N.O  flame  data.  The 
results  of  the  modeling  for  the  most  important  cases  studied  are  given  below  lor  the  flames  with  N  ()  as 
oxidizer. 


CQ-N^O  Flame  Structure 

Low  pressure,  laminar,  premixed  flames  of  C0-N20  have  been  stabilized  over  a  rectangular  flat-llame 
burner  (ref.  8).  Laser-induced  fluorescence  spectroscopy  was  used  to  establish  the  absence  of  CN,  CH, 
NH,  NH2  and  OH  in  these  flames.  Gas  chromatographic  sample  analysis  was  used  to  determine  the  CO, 
CO:,  N20,  NO,  N2  and  O,  concentration  profiles  for  three  C0-N20  flames  having  equivalent  ratios  of  1.00, 
1.32,  1.50.  Lean  flames  could  not  be  stabilized.  Temperature  profiles  for  all  three  flames  were  measured 
using  R  type  thermocouples.  Measured  temperature  profiles  were  corrected  for  radiation  losses.  These 
flames  are  considerably  lifted  above  the  burner  and  contain  a  single  luminous  zone.  The  main  feature  ot 
these  flames  is  the  absence  of  any  reactive  intermediates  except  oxygen  atom.  We  have  eliminated 
reactions  involving  the  element  H  and  reduced  the  272  step  mechanism  to  27  reactions.  This  mechanism 
was  then  used  to  model  these  C0-N20  flame  data. 

A  combined  "elementary  reaction  contribution"  and  "sensitivity"  analysis  showed  that  only  four  ol  the  27 
elementary  reactions  used  in  the  original  kinetic  mechanism  were  important  in  the  production  or  the 
consumption  of  major  species  in  the  flame.  A  comparison  of  the  calculated  concentration  profiles  using 
the  original  27-reaction  and  the  four-reaction  mechanisms  for  all  three  CO-N:0  flames  studied  in  this 
investigation  showed  that  the  difference  between  the  two  profiles  was  always  less  than  \%  of  the  total 
concentration.  Hence  we  propose  the  following  four-step  mechanism  for  the  CO-N-O  flames: 

CO  +  N20  —  CO;  +  N2 
n2o  +  M  -*  N2  +  O  +  *M 
nJo  +  O  —  N2  +  O, 

N20  +  O  —  NO  +  NO 

The  first  of  the  four  reactions  listed  above  is  the  most  important  reaction  for  C0-N,0  flames.  It  accounts 
for  almost  all  the  CO  consumption  and  nearly  90%  of  the  N;0  consumption.  The  second  reaction 
(N20+M— »N:+0+M)  is  an  important  reaction  for  kineticists.  It  plays  a  key  role  in  various  environments 
for  NO,  formation  or  decomposition.  It  is  also  used  for  generating  O  atoms  in  studies  of  elementary 
oxidation  reactions.  The  last  two  reactions  of  the  four-step  mechanism  listed  above  control  the 
concentration  profiles  of  02  and  NO. 


CH.-N-.Q  Flame  Structure 

Laminar,  premixed  flat  flames  of  CH4  with  N20  have  been  stabilized  and  studied  at  50  lorr  (ref.  5).  This 
study  represents  the  first  nearly  complete  study  of  the  structure  and  kinetics  of  CH4-N20  flames  including 
stable  and  unstable  species  measurements  and  detailed  chemical  kinetic  modeling.  Three  flames  were 
investigated,  with  slightly  fuel  rich,  near  stoichiometric  and  lean  mixtures.  Flame  modeling  has  been  done 
using  the  272  step  reaction  mechanism  and  the  reaction  mechanism  evaluated. 


The  results  of  the  flame  calculations  and  the  experiments  are  given  in  Figure  1.  The  comparison  shows 
good  quantitative  and  qualitative  agreement  between  the  measured  and  calculated  profiles  especially  for 
stable  species.  The  modeling  simulation  clearly  predicts  the  general  flame  structure  and  species 
concentration  profiles.  For  the  (lame  radicals,  the  theoretical  results  are  in  good  agreement  with  the 
experimental  results  for  NH,  CN  and  OH.  The  computed  maximum  for  the  CH  radical  is  shilled  in  the 
Hame.  The  exact  prediction  of  the  radical  concentrations  using  flame  modeling  is  difficult  owing  to  the 
high  reactivity  of  these  radicals,  especially  CH. 

The  reaction  mechanism  starts  with  the  thermal  decomposition  of  N,0  to  N,  and  O  which  is  the  initial 
reaction  for  radical  formation  in  the  flame.  The  oxygen  atom  formed  from  this  reaction  reacts  with  II, 
to  form  a  pool  of  H  and  OH  radicals  which  then  propagate  the  rest  of  the  mechanism.  The  results  of  the 
sensitivity  analysis  also  shows  that  N,0  is  not  totally  decomposed  to  N2,  but  it  forms  some  NH  and  NO 
molecules  directly  through  reaction  with  H  atoms.  The  mechanism  also  shows  that  methane  has  a  long 
induction  period  before  being  consumed  mainly  through  radical  reactions  (especially  with  OH  and  H)  to 
form  methyl  (CH,)  intermediate.  Once  methyl  intermediate  is  formed  it  has  two  reaction  paths  to  follow. 
The  first  path  is  to  form  formaldehyde  (CH:0)  which  then  starts  a  reaction  path  to  form  CO;.  The  ma  jor 
reaction  for  CO  conversion  to  C02  in  this  path  is  through  its  reaction  with  OH  radical.  The  second  path 
that  CH,  follows  is  to  form  CH2*  intermediate  which  then  forms  CH  through  radical  reactions,  The  CH 
formed  through  this  path  reacts  with  NO  from  N20  decomposition  to  produce  NH  or  CN  radicals.  This 
reaction  forms  the  first  link  in  hydrocarbon-nitrogen  interactions.  The  high  CN  and  NH  concentration  in 
the  flame  is  attributed  to  this  reaction.  However,  the  reaction  of  nitrous  oxide  (N20)  with  H  atoms  is 
another  important  source  for  the  formation  of  these  two  radicals  (CN  and  NH).  The  final  products  in  the 
exhaust  gas  were  mainly  N:,  H20,  CO  and  C02  with  small  amounts  of  NO. 


A.2.2  FLAMES  SUPPORTED  BY  NO,  AS  OXIDIZER 

Five  flames  supported  by  NO,  as  the  oxidizer  have  been  studied  recently  and  were  tested  against  our  272 
step  reaction  mechanism.  Volponi  and  Branch  (ref.  3  and  9)  studied  H,  and  C,H,  reacting  with  NO,  in 
an  argon  diluent.  They  measured  stable  species  concentrations  using  probe  sampling,  radical  species 
concentrations  using  LIF,  and  flame  temperatures  using  thermocouples  and  L1F.  Branch  et  al.  (rel.  ID) 
studied  two  different  flames  supported  by  NO,,  CHj/NOVO,  and  CH,0/N0:/0,.  They  measured  stable- 
species  concentrations  using  probe  sampling,  unstable  species  concentrations  using  LIF,  and  flame 
temperatures  with  a  thermocouple.  Zabamick  (ref.  11)  also  studied  the  CHj/N0,/0,  flame,  using  LIF  to 
measure  flame  temperatures  and  stable  and  radical  species  concentrations.  The  results  of  our  current 
modeling  of  the  hydrogen  and  the  acetylene  flames  are  most  indicative  of  the  reaction  mechanisms 
involved  and  are  discussed  below. 


H,-NQ,  Flame  Structure 

Measurements  of  the  composition  of  stable  and  unstable  species  and  temperature  in  laminar,  premixed, 
flat  flames  of  H,-NO,-Ar  have  been  made  and  compared  to  the  structure  calculated  with  a  flame  code 
including  detailed  chemical  kinetics  (ref.  9).  No  previous  detailed  flame  structure  measurements  and 
chemical  kinetic  modeling  of  this  flame  have  been  presented.  Similar  measurements  and  calculations  are 
reported  for  a  companion  H,-0,-Ar  flame  in  order  to  provide  a  comparison  to  a  previously  well 
characterized  flame.  We  have  modeled  these  flame  measurements  using  the  272  step  reaction  mechanism 
with  the  carbon  species  removed.  The  resulting  N7  step  mechanism  gives  good  agreement  to  the 
experimental  data  as  shown  in  Figure  2. 
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In  contrast  to  the  H:0-0:-Ar  name,  the  rate  and  sensitivity  calculations  for  H;,  NO,,  H,0  and  OH  in  the 
H,-NO,-Ar  flame  show  that  species  net  reaction  rates  are  usually  dominated  by  a  single  reaction.  Almost 
all  of  the  H,  consumption  and  H,0  lormation  is  by  0H+H,=H,0+H.  The  consumption  of  NO-,  and 
formation  of  NO  is  almost  entirely  by  NO,+H=NO+OH.  The  sum  of  these  two  reactions  gives  the  global 
reaction  of  the  name  H,+NO,=NO+H,0.  There  is  some  formation  of  O,  early  in  the  flame  by 
NO,+0=NO+0,  followed  by  O,  consumption  by  0,+H=0+OH.  The  OH  is  a  balance  of  formation  by 
NO,+H=NO+OH  and  consumption  by  0H+H,=H,0+H. 

The  reactions  with  greatest  negative  sensitivity  for  H,  and  NO,  have  the  highest  positive  sensitivity  tor 
NO  and  H,0.  The  reactions  with  the  largest  positive  sensitivity  for  H,  and  NO,  likewise  have  the  largest 
negative  sensitivity  for  NO  and  H,0.  The  OH  is  primarily  sensitive  to  its  major  formation  and 
consumption  reactions. 

The  importance  of  the  reaction  H,+NO,=HNO,+H  to  the  H,-NO,-Ar  flame  mechanism  is  observed  in  the 
sensitivity  and  analysis  for  H,,  NO,,  NO,  H,0  and  OH.  This  reaction  provides  the  most  important 
initiation  step.  It  also  has  a  significant  effect  on  the  OH  profile  by  producing  H  early  in  the  name  and 
promoting  the  formation  of  OH  by  NO,+H=NO+OH.  The  HNO,  is  subsequently  consumed  by  several 
reactions.  Using  a  mechanism  that  contains  only  the  three  reactions 

NO,  +  H  =  NO  +  OH 
OH'+  H,  =  H,0  +  H 
H,+  NO,‘  =  HNO,  +  H 

gives  a  flame  model  that  compares  to  within  5%  of  the  model  using  the  entire  reaction  mechanism. 


C,H,-NO,  Flame  Structure 

Volponi  and  Branch  (ref.  9)  studied  C,H,  reacting  with  NO,  in  an  argon  diluent  at  a  pressure  of  25  lorr. 
We  modeled  this  flame  using  the  272  step  reaction  mechanism  and  the  measured  temperature  profile.  The 
modeling  results  using  this  mechanism  are  close  to  those  Volponi  and  Branch  obtained  using  the  full.  331- 
reaction  mechanism  from  which  the  272-reaction  mechanism  was  derived.  This  indicates  that  modeling 
results  are  not  significantly  affected  using  the  smaller  mechanism. 

Measured  and  modeled  mole  fraction  profiles  for  major  species  are  shown  in  Figure  3.  Discrepancies  can 
be  seen  near  the  burner  surface,  as  the  modeling  overpredicts  the  surface  mole  fractions  of  unreaclcd 
species  and  underpredicts  the  surface  mole  fraction  of  product  species.  Agreement  is  good,  however,  in 
the  general  trends  shown  in  the  profiles,  and  in  the  final  mole  fractions  obtained  in  the  flame.  The 
modeling  was  used  to  find  the  principal  reactions  involved  in  each  specie's  production/consumption.  The 
principal  reaction  consuming  C,H,  is  C,H,+0H=C,H+H,0;  the  principal  reaction  consuming  NO,  is 
NO,+H=NO+OH;  the  principal  reaction  producing  NO  is  NO,+H=NO+OH;  the  principal  reaction 
producing  02  is  2N0,=2N0+0,;  the  principal  reaction  producing  CO,  is  NO,+CO=NO+CO and  the 
principal  reaction  producing  H,0  is  0H+C,H,=C,H+H,0.  The  importance  of  radicals  in  breaking  down 
the  unreactcd  species  is  evident. 

Comparisons  were  also  made  between  modeled  and  measured  mole  fractions  lor  the  radical  species  OH 
and  CN.  There  were  discrepancies  in  the  calculated  and  measured  trends  lor  these  species.  The  principal 
reaction  involving  OH  production  is  NO,+N=NO+OH,  and  the  principal  reaction  involving  CN  production 
is  HCN+0H=CN+H,0,  and  the  principal  reaction  involving  CN  production  is  HCN+0H=CN+H:0. 
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Volponi  and  Branch  (ref.  3)  also  modeled  C,H:  reacting  with  0:  with  Argon  gas  as  a  diluent  at  a  pressure 
ot  25  torr.  We  modeled  this  flame  using  the  272  step  mechanism  and  the  measured  temperature  profile. 
The  modeling’s  accuracy  was  on  a  par  with  what  was  obtained  for  the  CH./N-0  flame.  Modeling  was 
generally  in  good  agreement  with  experimental  data  for  major  species  concentrations  in  the  flames  with 
the  exception  of  H20,  whose  concentration  was  slightly  underpredicted.  Modeling  and  experiments  agreed 
in  the  general  trends  for  the  concentration  of  radical  species  OH  and  CH,  but  modeling  overpredicted  the 
peak  concentrations  of  both  species  by  a  factor  of  two.  As  with  the  CH./NO,  flame,  our  modeling  results 
compared  well  with  those  Volponi  and  Branch  obtained  using  the  331-reaction  mechanism. 


A.2.3  FLAMES  SUPPORTED  BY  NO  AS  OXIDIZER 

Zabarnick  (ref.  6)  studied  a  CH.,/N0/03  flame  at  63  torr,  using  LIF  to  measure  temperature  and  stable  and 
radical  species  concentrations.  This  flame  was  modeled  using  the  measured  temperature  profile  and  the 
272  step  mechanism. 

The  measured  NO  concentration  profile  showed  a  more  pronounced  drop-off  to  a  steady-stale  value  than 
modeling  predicted.  Also,  the  final  calculated  value  of  NO  in  the  flame  is  50%  higher  than  what  was 
measured;  NO  was  the  only  stable  species  Zabarnick  performed  measurements  on  in  this  work. 

Modeling  and  experimental  results  were  compared  for  radical  species.  The  predicted  peak  concentration 
of  CN  leads  the  experimentally  measured  peak  by  about  1  mm.  The  predicted  OH  concentration  showed 
a  drop-off  late  in  the  flame  that  was  not  detected  experimentally.  The  predicted  NH  peak  led  the 
measured  peak  by  about  .5  mm.  The  predicted  CH  peak  led  the  measured  peak  by  about  1  mm  also. 
Note  that  modeling  now  leads  the  experimental  data  for  the  location  of  peak  radical  concentrations,  while 
for  the  CR^N^O  flame  also  measured  by  Zabarnick  (ref.  6)  modeling  predictions  generally  lagged  the 
experimental  data. 


A.2.4  BURNING  VELOCITY  OF  FUEUNOx  LAMINAR  PREMIXED  FLAMES 

The  272  step  reaction  mechanism  was  used  to  model  the  burning  of  free-standing  flames  supported  by 
N:0  and  NO2,  and  the  modeled  flame  speed  was  compared  to  the  experimental  results  obtained  by  Parker 
and  Wolfhard  (ref.  12)  using  the  bunsen  burner  technique  (Table  2).  Also  shown  are  the  burning 
velocities  obtained  with  the  Miller  and  Bowman  (ref.  13)  mechanism  without  the  modifications  discussed 
in  the  present  investigation.  The  results  indicate  that  generally  good  agreement  is  obtained  between  the 
modeling  and  the  experimental  data  and  that  the  272  step  mechanism  gives  better  results  than  earlier 
mechanisms.  The  burning  velocity  of  the  flames  with  N:0  are  significantly  higher  than  those  with  NCL 
as  oxidizer.  Two  of  the  calculations  are  discussed  below  in  more  detail  in  order  to  indicate  the  importance 
of  fuel  and  oxidizer  decomposition  reactions  in  flame  propagation. 

Table  1  gives  the  calculated  and  measured  results  for  the  burning  rate  of  a  CH yN:0  flame.  The  calculated 
flame  velocity  is  90  cm/sec,  very  close  to  the  measured  velocity  of  105  cm/sec.  In  addition,  for  this  flame 
we  were  able  to  perform  analysis  on  what  reactions  contribute  most  strongly  to  a  given  specie's 
production/consumption.  It  was  found  that  the  major  reactions  leading  to  the  decomposition  ol  CH4  are 
the  following,  in  decreasing  order  of  their  rate  of  CH4  consumption: 


CH4  +  H  =  CHj  +  Hj 
CH4  +  OH  =  CHj  +  H,0 
CH4  +  O  =  CH,  +  OH' 

The  major  reactions  leading  to  the  breaking  down  of  N,0  were  found  to  be,  in  decreasing  order  of  the 
amount  of  N,0  they  consume: 

N,0  +  H  =  H,  +  OH 
n]o  +  M  =  N20  +  M 

The  importance  of  the  radicals  H,  OH  and  O  are  clearly  evident  in  the  propagation  of  this  (lame. 

A  free-standing  CjH4  flame  at  a  pressure  of  0.1  atm  was  modeled  and  the  calculated  flame  speed  was 
compared  to  that  measured  by  Parker  and  Wolfhard  (ref.  10).  The  results  are  included  in  Table  2,  and 
as  can  be  seen,  the  agreement  is  excellent  for  this  flame.  It  was  difficult  to  obtain  a  converged  solution 
modeling  many  Nonsupported  flames,  so  there  are  not  more  cases  reported  at  this  time.  It  is  possible 
that  the  larger  mechanism  used  to  model  this  flame  aided  in  obtaining  better  agreement  with  experimental 
data  than  was  obtained  in  several  of  the  cases  of  N,0-supported  flames.  More  cases  would  need  to  Ik 
run,  however,  to  draw  conclusions  on  the  relative  ability  of  the  two  mechanisms  to  model  free-standing 
flames  supported  by  NO,  and  N,0. 

Modeling  results  for  this  flame  showed  that  the  top  four  reactions  contributing  to  the  decomposition  of 
C,H4  in  the  flame  were  the  following: 

C,H4  +  OH  =  C,H,  +  H,0 
CH4  +  M  =  C,H,  +  H,  +  M 
C,[h4  +  H  =  CH,  +  H, 

C,H4  +  O  =  CH,  +  HCO 

The  top  four  reactions  leading  to  the  decomposition  of  NO,  were  found  to  be: 

NO,  +  N  =  NO  +  OH 
NO]  +  CO  =  NO  +  CO, 

NCt]  +  M  =  NO  +  0  +  *M 
No]+  OH  =  HO,  +  NO 

Initiation  of  the  chemistry  in  these  flames  can  be  seen  to  be  closely  tied  to  reactions  producing  radical 
species. 


A.2.5  SUMMARY  AND  CONCLUSIONS 

We  have  completed  detailed  comparisons  of  calculations  and  measurements  of  the  structure  and  burning 
velocity  of  fucl/NO,  mixtures  using  a  272  step  reaction  mechanism.  The  comparison  between  the 
calculated  flame  structure  and  the  experimental  flame  structure  for  stable  species  was  found  to  be  very 
good  for  fucl/N,0  flames  and  good  for  fuel/NO,  flames.  The  concentration  profiles  for  radical  species 
were  found  to  be  generally  well  represented  qualitatively  but  not  well  represented  quantitatively.  It  was 
concluded  that,  despite  some  remaining  difficulties  with  the  reaction  mechanism,  it  appeared  to  be  reliable 
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in  describing  the  overall  combustion  behavior  of  a  wide  range  of  fuel  and  oxidizer  mixtures. 

The  most  important  reactions  of  the  oxidizer  (N:0  or  NO:)  are  with  H  atoms  and,  to  a  lesser  extent,  with 
CO.  Reaction  of  either  oxidizer  with  H  is  a  chain  propagating  reaction,  in  contrast  with  the  chain 
branching  reaction  of  H  with  0;  which  is  of  equal  importance  in  fuel  oxidation  by  O,.  In  addition,  the 
reaction  of  NO,  with  H  is  slower  than  N,0  with  H.  Finally,  the  reaction  of  N;0  with  CO  can  be 
significant  both  in  consumption  of  CO  and  formation  of  CO,.  This  situation  is  again  in  contrast  with  the 
oxidation  chemistry  of  systems  by  O,  in  which  the  conversion  of  CO  to  CO,  is  almost  entirely  be  reaction 
of  CO  with  OH.  The  difference  between  the  use  of  N,0  or  NO,  as  oxidizer  is  that  the  former  produces 
primarily  N,  while  the  latter  produces  NO.  The  subsequent  slow  reduction  of  NO  to  N,,  even  when  it  is 
thermodynamically  favored,  accounts  for  the  most  striking  difference  between  the  two  oxidizers.  The  most 
important  effects  of  the  oxidizer,  therefore,  are  that  the  nitrogen  oxides  are  less  effective  chain  carriers 
and  lead  to  slower  reaction  rates  compared  to  O,. 

The  oxidation  reactions  of  CH4  in  the  presence  of  NOx  is  generally  similar  to  the  oxidation  by  O,.  Chain 
propagating  and  branching  reactions  lead  to  the  formation  of  H,  O  and  OH  and  these  species  progressively 
abstract  hydrogen  and  partially  oxidize  CH3  to  CH,0.  The  CH,0  is  then  converted  largely  to  CO  through 
the  intermediate  HCO.  Subsequent  oxidation  of  CO  to  CO,  is  by  reaction  with  NOx  of  OH  as  mentioned 
above.  If  the  fuel  is  CH,0  instead  of  CH4,  the  latter  stages  of  this  chain  become  dominant.  If  the  fuel 
is  C,H,  instead  of  CH4,  then  reactions  of  CH,  become  more  important  and  CH,0  becomes  less  important. 
Rate  constants  for  the  reactions  involved  in  the  hydrocarbon  chemistry  derived  from  previous  studies  are 
generally  successful  in  the  description  of  the  transformation  of  CH,  to  CO. 

The  third  major  aspect  of  fuel/NOx  flame  chemistry  is  the  interaction  of  H,  C  and  N  containing  species. 
Reaction  of  CH,  species  with  NO  to  form  HCN  and  subsequent  reaction  of  HCN,  CN  or  NCO  with  NO 
lead  to  the  conversion  of  NO  to  N,.  This  process  is  obviously  of  most  importance  when  NO,  or  NO  is 
the  oxidizer  rather  than  when  N,0  is  oxidizer,  since  in  the  latter  case  N,  is  formed  directly.  This  scheme 
also  shows  the  essential  features  of  the  oxidation  process  in  the  case  where  HCN  is  formed  as  a  lucl 
during  the  decomposition  of  the  energetic  solid. 
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TABLE  1.  Rate  constants  used  in  the  present  272  Step 


Reaction  Mechanism  which  differ  from  those  of 
Volponi  and  Branch  (Ref.  3).  Units  are  mole,  cm’,  sec, 
K,  and  Cal/mole. 

Reaction 

A 

n 

Ea 

N.O+OH«Nj+HOj 

1.00el3 

0 

10000 

N;0+H=N;+OH 

1.50el4 

0 

15090 

NjO+M+Nj+O+M 

4.90elS 

0 

57500 

N.0+0*Nj+05 

7.00el4 

0 

28200 

N.0+0*2N0 

5.60el4 

0 

28200 

CO+N-OaCO.+NO 

2.00c  12 

0 

17300 

TABLE  —  Comparison  of  measured  and  calculated 
flame  speed  of  fuel/NO,  mixtures*. 


REACTANTS  MEASURED*  CALCULATED 

M/B*  This  Wort* 


Hj/N-0 

300 

227 

240 

OL/N;0 

105 

162. 

90 

C,IVN;0 

160 

128 

150 

C,H./N;Q 

110 

94 

102 

C-HVNOj 

623 

61 

61 

(a)  All  flame  speeds  given  in  cm/sec  at  0.10  atm  and  for 
a  stoichiometric  mixture. 

(b)  Parlcer  it  Wolfhard  (Ref.  12). 

(c)  Calculation  using  mechanism  of  Miller  &  Bowman 

(Ret.  13). 

(d)  Calculation  using  the  present  272  step  mechanism. 


Flam*  petition  (mm) 

Figure  1.  Comparison  of  CH*/NjO  flame  structure  data 
of  Habeebullah  el  al.  (Ref.  5)  to  calculations  of  flame 
structure  using  272  step  reaction  mechanism. 


Tlinplialuii  (K) 


Position  aoovs  mo  oumor  (cm) 


Figure  2.  Comparison  of  HyNQ,  flame  structure  data 
of  Volponi  and  Branch  fRef.  9)  to  calculations  of  Flame 
Structure  Using  772  step  reaction  mechanism. 


Figure  3.  Comparison  of  C,iyNO,  flame  structure 
data  of  Volponi  and  Branch  (Ref  3)  to  calculations  of 
Dame  structure  using  772  step  reacuoe  mechanism. 


B.  FUNDAMENTAL  STUDIES  OF  ROCKET  COMBUSTION  INSTABILITY 
B.l  RESEARCH  OBJECTIVES 


Rocket  engine  manufacturers  rely  heavily  on  empirical  and  actual  test  data  during  the  design  stage  of 
rocket  engines  in  order  to  prevent  the  occurrence  of  instabilities  during  (light  operation.  This  is  clearly 
an  expensive  approach.  Thus,  the  long  term  objective  of  our  research  is  to  develop  predictive  models  that 
will  help  rocket  designers  build  high  performance,  stable  liquid  rocket  engines. 

A  complex  interaction  of  several  processes  such  as  atomization  of  the  primary  fuel  jet,  secondary 
atomization  of  the  resulting  drops,  vaporization,  mixing  and  combustion  are  responsible  for  liquid  rocket 
engine  instability.  Since  atomization  is  the  first  mechanism  in  the  chain  of  events  leading  to  eventual 
combustion  and  possibly  combustion  instability,  our  near  term  objective  is  to  understand  the  influence 
of  atomization  on  the  stability  properties  of  rocket  engines.  The  focus  of  our  work  is  to  elucidate  and 
model  the  response  of  atomization  under  controlled  modulated  conditions. 


B.2  BACKGROUND 

Combustion  instability  that  is  characterized  by  large  amplitude  combustion  chamber  pressure  oscillations 
is  a  commonly  observed  feature  in  many  liquid  propelled  rocket  engines.  Severe  engine  damage  could 
result  due  to  mechanical  failure  and/or  enhanced  heat  transfer  to  the  chamber.  Liquid  fuel  injected  into 
a  rocket  engine  combustion  chamber  undergoes  a  sequence  of  processes  prior  to  exiting  the  chamber  in 
the  form  of  products  of  combustion.  A  typical  scenario  leading  to  combustion  instability  is  as  follows: 
an  initial  chamber  pressure  disturbance  induces  fluctuations  in  atomization  rale.  This  results  in  variation 
of  propellant  droplet  size  distribution,  which  affects  vaporization,  mixing  and  reaction  rates,  inducing 
fluctuations  in  the  engine  burning  rate.  This  process  can  feed  energy  back  into  the  pressure  field,  thus 
closing  the  loop.  If  proper  phase  relationship  is  achieved  during  these  processes,  amplification  of  initial 
disturbances  can  occur,  leading  to  combustion  instability.  Integration  of  all  the  physical  processes 
involved  into  a  comprehensive  model  of  engine  stability  is  not  probable,  unless  each  process  is  well 
understood. 

Injectors  in  liquid  rocket  engines  are  exposed  to  relatively  high  amplitude  pressure  and  velocity  field 
oscillations  (acoustic  waves).  These  oscillations  have  a  pronounced  effect  on  the  atomization  of  liquid 
propellant.  Sufficient  evidence  from  past  work  and  our  present  research  reveals  that  the  structure  and 
characteristics  of  atomization  under  unsteady  conditions  (as  would  occur  during  a  combustion  instability) 
deviate  considerably  from  atomization  behavior  under  steady  conditions.  Thus,  the  dynamic  behavior  of 
liquid  propellant  atomization  is  of  major  importance,  especially  when  atomization  produces  very  small 
droplets. 

Miessc  (1955)  showed  that  significant  acoustic  interactions  can  be  achieved  between  a  low  speed  jet  and 
local,  high  amplitude  sound  fields.  Reba  and  Brosilow  (1960)  conducted  the  first  experiment  on  liquid 
jet  behavior  perturbed  by  a  longitudinal  acoustic  wave  created  by  a  siren  (their  work  is  also  reported  by 
Torda  and  Schmidt,  1964).  Buffum  and  Williams  (1967)  investigated  the  behavior  of  a  liquid  jet  subject 
to  a  transverse  standing  acoustic  field.  They  studied  the  periodic  oscillations  of  the  liquid  jet  position  in 
space  and  argued  that  this  phenomenon  might  explain  combustion  instability  in  liquid  rocket  engines  with 
impinging  jet  type  injectors.  Jet  oscillatory  deviations  cause  fluctuations  in  impingement  point  and 
geometrical  injection  parameters,  which  can  induce  variations  in  the  atomization  rate  and  pattern.  Their 
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results  showed  a  strong  correlation  between  the  acoustic  field  and  the  jet  behavior  for  forcing  amplitudes 
of  160  dB  and  frequencies  up  to  500  Hz. 

Heidmann  (1965)  studied  acoustic  field  and  jet  interaction  in  a  two-dimensional  circular  combustor.  He 
found  a  strong  correlation  between  the  acoustic  field  and  the  oscillatory  jet  breakup  behavior,  Ingebo 
(1966)  studied  the  behavior  of  liquid  jet  atomization  perturbed  by  longitudinal  acoustic  waves.  A  siren 
periodically  obstructed  the  chamber  exit  and  created  longitudinal  acoustic  perturbations.  Streak 
photography  showed  reversal  of  the  gas  flow  during  instability.  High  speed  photography  provided  the  first 
measurement  of  drop  size  distribution  under  unsteady  and  hot  environment.  He  showed  that  under 
resonant  conditions,  atomization  and  vaporization  rates  are  significantly  enhanced.  The  jet  breakup  length 
was  one  sixth  compared  with  the  unforced  case.  Much  finer  droplet  size  was  obtained  (62  microns  in 
resonant  forcing  conditions  and  217  microns  for  unforced  cases).  The  length  for  complete  vaporization 
decreased  from  -50  cm  to  5  cm.  More  recently,  Lecourt  and  Foucaud  (19S7)  studied  the  stability  of 
several  injectors  using  the  forcing  method  developed  by  Ingebo.  Due  to  extreme  operating  conditions 
(1000  psia),  only  pressure-time  traces  were  recorded. 

Heidmann  and  Groeneweg  (1969)  studied  the  transfer  function  between  pressure  oscillations  and 
atomization  rate.  They  argued  that  if  the  in-phase  transfer  function  is  positive,  unstable  combustion  might 
occur.  The  justification  for  this  approach  comes  from  the  Rayleigh  criterion.  Different  types  of  transverse 
acoustic  oscillations  (travelling,  standing,  radial  transverse  mode  and  vortex  mode,  characterized  by  a 
steady  tangential  velocity  component)  were  investigated  experimentally  and  analytically.  They  found  that 
variations  in  the  peak  values  of  the  response  factors  are  caused  by  the  velocity  sensitivity  of  the 
atomization  process,  and  a  velocity  perturbation,  in-phase  with  pressure,  elevated  the  peak  value  and  is 
destabilizing  to  a  combustion  system.  However,  an  out-of-phase  velocity  perturbation  reduced  the  peak 
value  of  the  response  factor  and  is  stabilizing.  They  predicted  that  potentially  the  maximum  average 
response  for  the  traveling  mode  is  larger  than  for  either  the  standing  or  radial  modes.  Both  their  analytical 
and  experimental  results  indicated  that  a  steady  vortex  velocity  had  an  amplifying  effect  and  increased  the 
peak  value  of  the  response  factor  by  an  order  of  magnitude. 

Overall,  very  little  is  known  about  the  dynamics  of  atomization  under  modulated  conditions.  Current 
interest  in  combustion  instabilities  of  liquid  rocket  engines  and  the  advent  of  modern  diagnostics  has 
initiated  our  study  of  atomization  behavior  under  oscillatory  environment.  In  our  work  external  forcing 
is  used  to  simulate  combustion  induced  pressure  and  velocity  fluctuations.  The  advantage  is  a  complete 
control  over  frequency,  and  amplitude,  so  that  linear  and  nonlinear  forcing  regimes  can  be  investigated. 
Our  study  focuses  on  the  phase  and  amplitude  relationship  between  forcing  and  atomization  response. 
Although  these  quantities  have  been  shown  to  be  important  in  theoretical  studies  on  the  possible  coupling 
between  atomization  and  acoustic  waves  leading  to  combustion  instability,  they  have  never  been  measured. 
This  knowledge  can  lead  to  the  use  of  linear  stability  theory  to  clarify  the  role  of  atomization  in  high 
frequency  combustion  instability. 


B.3  DESCRIPTION  OF  EXPERIMENTAL  FACILITY 

As  one  of  the  objectives  of  this  research  program  is  to  provide  useful  atomization  information  tor  the 
design  of  actual  injectors,  it  is  important  to  give  special  attention  to  similarity  parameters.  Moreover,  it 
is  extremely  difficult  to  model  a  real  rocket  engine  because  of  the  severe  conditions  encountered  in  an 
.actual  combustion  chamber.  High  pressure  and  intense  heal  transfer  rates  precludes  accurate  measurement 
of  quantities  such  as  pressure,  temperature  and  droplet  size  distribution.  Propellants  used  in  the  rocket 
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industry  are  usually  toxic  and  hypergolic.  The  approach  adopted  is  to  identity  dimensionless  numbers, 
which  govern  atomization  and  try  to  maintain  them,  or  at  least  operate  in  the  proper  ramie.  Four 
dimensionless  numbers,  liquid  jet  Reynolds  number 


Weber  number  based  on  gas  phase  properties 


/V/V  and  P/P,,  are  used  to  characterize  atomization!  Here,  the  symbols  /<  and  p  represent  dynamic 
viscosity  and  density.  The  subscripts  l  and  g  indicate  liquid  and  gas.  However,  the  definition  of  the 
quantities  which  are  chosen  to  calculate  the  Reynolds  and  Weber  numbers  varies  in  the  literature. 

The  turbulence  level  in  the  liquid  phase,  ns  the  jet  emerges  from  the  nozzle  is  important  because  it 
characterizes  the  amplitude  of  the  disturbances  on  the  jet  surface,  which  are  likely  to  be  amplified  by 
dynamic  interactions  with  the  surrounding  gas.  This  justifies  our  definition  of  the  Reynolds  number.  The 
Weber  number  based  on  gas  phase  properties  takes  into  account  aerodynamic  forces,  the  main  mechanism 
for  jet  breakup  and  atomization  at  high  velocities.  Although  no  single  mechanism  can  explain  the  entire- 
process  of  atomization,  aerodynamic  interactions  between  the  liquid  jet  and  the  surrounding  gas  provide 
the  main  destabilizing  agent  (Reitz  and  Bracco,  1986). 

A  survey  of  studies  on  coaxial  injectors  has  been  completed.  The  operating  range  of  relevant  results  is 
reported  in  a  Re,  -  Weg  diagram  in  Figure  1.  The  two  horizontal  lines  correspond  to  theoretical  results 
from  Ranz  (1958)  and  Miesse  (1955),  who  predicted  that  atomization  regime  lies  beyond  a  critical  We.. 
Our  approach  is  to  operate  in  the  atomization  regime  and  come  as  close  as  possible  to  the  configuration 
under  simulation  by  a  proper  choice  of  simulants. 

B.3.1  Simulation  Selection 

Since  real  atomizer  injection  velocities  are  very  high,  they  cannot  be  increased  to  maintain  similarity.  A 
survey  of  different  propellant  combinations  revealed  that  the  best  which  can  be  done  is  to  maintain  the 
value  of  liquid  and  gas  velocity.  In  order  to  maintain  geometrical  similarity,  we  have  retained  the  inner 
post  diameter  (Figure  2).  Therefore,  the  dimensionless  numbers  indicate  that  the  simulant  selection  is 
governed  by  the  following  two  equations. 


(LI)  *  (Ll\ 

1  'actual  I  'model 


actual 


model 
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Tables  1  and  2  show  the  values  of  p Jft,  and  p Ja,  lor  different  combinations  for  both  hot  and  cold  run 
simulations  of LOJCH,  rocket  propellants.  For  cold  How  testing,  the  combination  Freon  113,  RI3B  seems 
adequate,  if  atomization  data  on  LOJCH :  is  to  be  obtained.  For  hot  How  tests,  one  sees  that  tor  all 
simulants  (hydrogen  peroxide,  ethanol)  the  value  of  falls  short  by  a  factor  of  at  least  2. 

In  our  case,  we  chose  to  use  a  combination  of  water  and  air  mixed  with  refrigerant  for  cold  flow  runs. 
For  experiments  with  combustion,  the  ethanol/air  mixed  with  oxygen  combination  has  been  selected. 

B.3.2  Combustor  Design 

The  size  of  the  combustor  chamber  has  to  be  chosen  carefully  to  ensure  almost  complete  combustion 
before  the  nozzle  exit.  In  particular,  considerations  using  parameters  such  as  L'  or  tu  helped  us  in  our 
task.  The  quantity  L"  is  a  ratio  of  chamber  volume  to  throat  area.  Values  of  L'  must  be  chosen  in 
accordance  with  the  fuel-oxidizer  choice.  The  fluid  residence  time  in  the  chamber  is  denoted  by  /i(.  Also, 
the  selected  excitation  system  places  a  bound  on  the  size  of  the  chamber,  because  large  amplitude 
modulations  are  possible  only  if  the  chamber  volume  is  small  enough. 

The  combustion  chamber  design  is  shown  in  Figure  3.  The  uni-element  coaxial  combustor  has  transverse 
dimensions  of  5  cm  x  5  cm.  The  chamber  is  designed  to  withstand  a  pressure  of  10  atm  absolute.  This 
limitation  comes  from  the  quartz  windows  (1.25  cm  thick).  The  chamber  length  can  be  varied  from  2<S 
to  43  cm  by  moving  the  exit  nozzle  location.  The  adjustable  length  allows  for  the  study  of  various 
combustion  efficiencies  and  longitudinal  acoustic  modes.  The  length  of  the  chamber  dictates  the  residence 
time  of  the  reactants  and  hence  the  amount  consumed,  or  in  essence  the  combustion  efficiency.  Two 
quartz  windows,  25  x  5  cm,  on  opposite  sides  of  the  chamber  provide  full  optical  access  from  the  injection 
plane  to  25  cm  downstream.  The  operating  pressure  of  the  chamber  for  given  fuel  and  oxidizer  mass  flow 
rates  can  be  varied  by  using  exit  nozzles  of  different  throat  area.  The  chamber  is  mounted  vertically  to 
eliminate  liquid  accumulation  during  cold  flow  tests.  The  injector  head  is  made  of  a  5  x  5  x  7  cm 
stainless  steel  block,  as  shown  in  Figure  2.  The  injector  head  can  easily  be  modified  to  accommodate 
other  types  of  atomizers  such  as  an  impinging  jet.  The  center  post  of  the  injector  is  machined  to  allow 
for  smooth  reduction  of  liquid  fluid  passage  area.  The  mounting  of  the  center  post  can  be  adjusted  to 
achieve  the  desired  recess  relative  to  the  injection  plane.  Center  posts  of  different  design  and  diameter 
can  be  mounted  into  the  injector  block.  The  injector  used  in  this  study  has  a  center  post  with  inside 
diameter  of  1  mm  with  a  4  mm  width  of  the  surrounding  annulus. 

B.3.3  Excitation  Method 

Acoustic  compression  drivers  are  used  to  excite  transverse  modes  in  the  combustion  chamber.  The 
experimental  configuration  is  shown  in  Figure  4.  This  kind  of  perturbation  is  chosen  because,  in  many 
high  frequency  combustion  instability  cases,  the  coupling  between  combustion  and  chamber  transverse 
acoustic  modes  have  been  recorded.  The  use  of  acoustic  excitation  devices  in  studies  in  the  field  of  active 
control  of  combustion  instability  has  proved  to  be  successful.  In  particular,  very  high  pressure  excitation 
levels  can  be  obtained. 


Most  acoustic  drivers  operate  adequately  in  the  frequency  range  from  1()(K)  Hz  to  boot)  Hz.  Moreover, 
many  high  frequency  combustion  instabilities  are  encountered  in  this  range.  Hence,  the  excitation 
frequency  domain  for  transverse  excitation  is  chosen  as  above.  This  sets  the  chamber  transverse 
dimension  to  under  5  cm.  It  is  also  desirable  to  excite  the  chamber  at  a  frequency  different  front  any  of 
its  characteristic  frequencies,  in  order  to  avoid  resonance.  Thus,  we  can  exercise  complete  control  on  the 
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amplitude  of  modulation.  High  forcing  pressure  amplitudes  can  be  obtained  by  tuning  the  length  of  the 
ducts  connecting  the  compression  chambers  to  the  combustion  chamber. 


B.4  DIAGNOSTIC  TECHNIQUES 

Diagnostic  techniques  aimed  at  studying  this  configuration  are  under  development.  In  particular,  the  time 
variation  of  quantities  such  as  atomization  rate,  jet  breakup  length,  length  for  complete  atomization,  under 
forcing  and  non-forcing  conditions  are  of  interest.  Large  and  small  scale  fluid  mechanical  fluctuations 
and  their  phase  with  respect  to  external  excitation  are  also  important  in  belter  understanding  the  dynamic 
behavior  of  atomization.  The  task  is  difficult  because  the  frequency  range  of  interest  is  above  1000  Hz. 
Currently,  phase  locked,  back-lighted  photography  and  Schlieren  arc  used.  Diagnostic  systems  like  the 
ones  in  Figures  5  and  6  can  give  access  to  the  aforementioned  interesting  quantities.  The  signal  from  the 
acoustic  drivers  is  used  to  periodically  trigger  a  light  source.  An  optical  arrangement  focuses  the  light 
onto  a  region  of  the  atomizing  jet.  The  resulting  image  is  captured  by  a  continuous  high  speed  camera. 
Phase  averaging  is  then  used  to  minimize  cycle-to-cycle  variations.  The  images  arc  then  processed  to 
extract  jet  breakup  length  and  characteristic  time  for  complete  atomization.  By  shilling  the  phase  of  the 
light  trigger  with  respect  to  the  excitation  signal,  phase  sensitive  images  of  the  atomizing  jet  arc  obtained. 

Currently,  the  response  of  the  chamber  pressure  and  localized  heat  release  in  the  combusting  region  to  an 
acoustic  perturbation  is  under  investigation.  In  particular,  the  amplitude  and  phase  of  the  response  relative 
to  the  perturbation  is  being  analyzed.  A  piezoelectric  transducer  mounted  on  the  chamber  wall  measures 
the  chamber  pressure  response.  This  type  of  transducer  is  capable  of  measuring  high  frequency  and  low 
amplitude  acoustic  waves.  The  localized  heat  release  is  measured  by  collecting  the  light  emitted  by  C2 
radicals,  which  are  only  present  in  the  high  temperature  combusting  regions  of  the  flow.  The 
measurement  is  local  in  the  sense  that  an  optical  system  only  gathers  light  from  a  very  small  region  in  the 
combustion  region.  The  emitted  light  is  collected  and  focused  onto  a  photomultiplier  lube  for  signal 
generation.  In  this  manner,  a  relative  measure  of  the  localized  heat  release  can  be  obtained  for  comparison 
between  excited  and  unexcited  conditions. 

In  addition,  laser  back  scattering  and  extinction  diagnostic  systems  are  under  development,  as  shown  in 
Figure  7.  These  systems  are  capable  of  measuring  the  atomization  response  to  an  imposed  perturbation. 
The  idea  is  to  focus  a  He-Ne  laser  beam  into  the  initial  atomizing  region  near  the  inner  post  lip,  the 
scattered  light  is  then  collected  by  a  series  of  PMT  detectors  and  the  extinction  of  the  laser  beam  is 
measured  by  a  photo  diode.  The  intensity  of  the  laser  beam  received  by  the  photo  diode  is  directly 
proportional  to  the  amount  of  fluid  which  has  passed  through  it.  Spectral  analysis  of  the  signals  can  then 
give  information  concerning  the  modification  of  frequency  content  of  the  atomization  process  due  to 
acoustic  perturbations. 


B.5  RESULTS 

B.5.1  Cold  Flow  Results 

A  scries  of  pressure  measurements  inside  the  chamber  have  been  conducted  in  order  to  obtain  the 
fluctuating  pressure  and  velocity  field  generated  by  the  excitation  system.  The  measurements  indicate  thni 
transverse  acoustic  field  strengths  of  up  to  156  dB,  and  acoustic  velocity  fluctuations  of  up  to  4.2  m/s  at 
the  injector  location  can  be  produced. 
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A  numerical  code  using  linear  acoustic  equations  has  been  developed  to  model  the  pressure  field  in  the 
test  section.  The  equations  are  solved  using  a  second  order  accurate  finite  difference  scheme.  Input 
includes  the  acoustic  driver  boundary  conditions  and  a  temperature  map  of  the  system  to  account  for  the 
speed  of  sound  distribution.  Output  includes  the  acoustic  modes  for  the  system,  and  resulting  pressure 
and  velocity  fields.  The  code  will  be  enhanced  to  include  flow  within  the  chamber  to  belter  model  real 
conditions. 

A  set  of  experiments  has  been  carried  out  to  determine  the  effect  of  acoustic  forcing  on  the  global 
behavior  of  the  fuel  spray.  These  measurements  are  difficult  due  to  the  accumulation  of  fluid  on  the 
windows  of  the  chamber.  The  results  indicate  that  under  out-of-phase  forcing  conditions,  there  is  an 
increase  in  the  volume  of  fuel  droplets  drawn  into  the  recirculating  region  at  the  top  of  the  chamber. 
Visually,  the  level  of  fluid  accumulated  on  the  window  rises  with  excitation.  Further  experiments  are 
planned  to  elucidate  the  physical  phenomenon  and  obtain  precise  quantitative  results. 

B.5.2  Hot  Flow  Results 

The  combustion/excitation  system  has  been  improved  and  made  ready  for  experimentation.  Smooth 
reliable  combustion,  as  well  as  integration  with  the  excitation  system  has  been  obtained.  The  operating 
limits  involving  reactant  parameters  and  chamber  geometrical  settings  have  been  explored.  In  addition, 
the  diagnostic  Sclierien  imaging  technique  has  been  refined  to  generate  high  qualily(eonlrast)  imaging  of 
the  flow  during  hot  runs,  as  shown  in  Figure  S. 

Recent  experimentation  has  been  conducted  to  determine  the  effect  of  equivalence  ratio  and  excitation 
amplitude  upon  the  time  variation  of  the  chamber  pressure  between  forced  and  unforced  conditions.  The 
gas  phase  and  liquid  phase  velocities  were  held  constant  at  92  m/s  and  15  m/s,  respectively.  The 
equivalence  ratio  was  varied  by  the  addition  of  pure  oxygen  to  the  gas  phase  resulting  in  a  range  from  3.S 
to  7.4  .  These  procedures  were  carried  out  in  order  to  keep  the  atomization  characteristics  of  the  system 
constant.  The  enrichment  of  the  gas  phase  with  oxygen  increased  the  gas  phase  velocity  very  little. 
Moreover,  the  aerodynamic  properties  of  oxygen  are  essentially  similar  to  air.  Thus,  atomization 
parameters  including  Weber  number  and  density  ratio  were  held  relatively  constant.  However,  the  static 
chamber  pressure  varied  with  equivalence  ratio.  This  was  caused  by  the  use  of  the  same  size  chamber 
exit  nozzle  for  all  the  experiments.  A  variety  of  exit  nozzles  are  under  construction  in  order  to  decouple 
the  chamber  pressure  from  equivalence  ratio.  As  an  example,  for  an  experiment  with  an  equivalence  ratio 
near  one,  which  is  a  high  heat  release  and  hence  high  chamber  pressure  case,  a  larger  exit  nozzle  will 
be  used  to  reduce  the  chamber  pressure  so  that  it  will  remain  constant. 

Unexcited  chamber  pressure  recordings  for  equivalence  ratios  of  7.4  and  3.8  are  shown  in  Figures  9  and 
13,  respectively.  The  power  spectra  of  these  signals,  shown  in  Figures  10  and  14,  indicate  the  presence 
of  three  frequency  ranges.  The  first  range  is  centered  at  200  Hz  and  the  other  two  ranges  are  centered 
at  1250  Hz  and  2500  Hz.  There  is  a  slight  increasing  shift  of  these  ranges  as  the  equivalence  ratio  is 
reduccd(approaches  one).  This  is  probably  caused  by  the  fact  that  as  stoichiometric  conditions  art- 
approached,  the  chamber  experiences  higher  temperatures  and  thus  higher  sound  speeds,  and  hence 
increased  oscillation  frequency. 

The  acoustic  compression  drivers  were  tuned  to  excite  the  chamber  at  295(1  Hz.  This  frequency  was 
chosen  because  it  corresponds  to  a  largo  amplitude  pressure  fluctuation  in  the  chamber  (determined 
experimentally),  and  it  is  in  the  high  frequency  (screeching)  range.  Chamber  pressure  records  under  forced 
conditions,  for  equivalence  ratios  of  7.4  and  3.8  are  shown  in  Figures  1 1  and  15,  respectively.  The  power 
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spectra  of  these  signals,  shown  in  Figures  12  and  16,  indicate  a  peak  at  the  excitation  frequency.  This 
is  expected  because  the  combustion  chamber,  excitation  system,  and  pressure  transducer  are  in  a  closed 
system.  The  interesting  result  is  that  the  observed  peak  at  the  excitation  frequency  increases  as  the 
equivalence  ratio  is  decreased.  The  peak  power  at  the  excitation  frequency  versus  equivalence  ratio  is 
displayed  in  Figure  17.  This  result  is  biased  by  the  fact  that  the  static  chamber  pressure  also  increased 
with  decreasing  equivalence  ratio.  As  stated  earlier,  the  elimination  of  this  bias  is  underway. 

The  effect  of  excitation  amplitude  for  a  fixed  equivalence  ratio  of  3.8,  holding  all  other  parameters 
constant,  was  examined.  As  shown  in  Figure  18,  the  peak  power  at  the  excitation  frequency  increases  as 
the  excitation  amplitude  is  increased.  Again,  this  was  expected  because  the  pressure  transducer  and 
excitation  system  are  in  a  closed  system. 

Currently,  work  is  under  way  to  determine  the  effect  of  gas  phase  velocity,  equivalence  ratio,  excitation 
frequency,  chamber  length,  and  injector  post  recess  upon  the  time  variation  ol  chamber  pressure  and 
localized  heat  release  (C2  emission)  between  forced  and  unforced  conditions.  Schlierien  images  are  being 
acquired  to  analyze  the  fluid  mechanical  response  to  excitation,  i.e  (jet  breakup  length,  length  for  complete 
atomization,  vortical  structure,...etc).  The  analysis  includes  an  examination  of  the  frequency  content 
(power  spectra)  of  the  measured  signals  to  determine  it  the  acoustic  lorcing  frequency  can  drive  the 
physical  process  under  investigation.  The  phase  and  amplitude  response  ol  the  signals  relative  to  the 
perturbation  will  also  be  obtained.  This  information  is  a  direct  measure  of  the  (n.T)  parameters  applied 
to  (n,x)  correlations  commonly  used  in  industry  to  predict  instability  behavior.  In  addition,  auto-  and 
cross-correlations  of  these  signals  will  be  acquired  to  determine  the  periodicity  and  similarity  of  the 
measured  signals  and  perturbation. 


B.6  FUTURE  WORK 

Laser  backscattering  and  extinction  experiments  are  planned  to  obtain  real  time  data  on  spray  fluctuations 
in  order  to  measure  the  frequency  response  of  atomization  to  acoustic  perturbation.  Furthermore,  a  Phase 
Doppler  system  is  under  development  in  order  to  measure  the  fluctuation  in  drop  size  and  drop  velocity 
in  the  atomizing  region.  Thus,  this  system  will  provide  precise  quantitative  measurements  of  real  physical 
atomization  parameters  and  give  a  more  clear  understanding  of  atomization  response. 

In  addition,  an  image  processing  system  has  been  acquired.  This  system  will  generate  sharp,  detailed 
images  of  the  spray  and  flow  structure.  Phase  averaging  can  then  provide  statistically  reliable 
quantitative  information  concerning  fluid  mechanical  response. 
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Table  1:  pi  l  pi  of  Liquids 


LOX  (crit) 

LOX  (sat) 

Water 

Ethanol  |  Freon  113 

Jet  A 

Shell  iv ax 

worn 

2.6 

7.2 

1.0 

1.1 

2.3 

0.5 

0  4 
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Table  2:  pj/crj  of  Different  Gas  Liquid  Combinations 


|  LOX(crit)/GH2  | 

LOX(sat)/GH2 

HiOil propane 

Ethanol/Oj 

Water/  R13B  ' 

m 

1176 

375  |  80  j 

195 

216  , 

i  Freonll3/R13b  |  HjO j/  propane.  R13B  |  ethanol/Oj,  Air  j  Water,  ethanol/  R.13B 


Figure  3.  Picture  u  CurtbusMon  Chamber 
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Figure  4:  Combustor  Sketch  for  Transverse  Mode  Forcing 


opocjl  tyncm  Gun  human  chamber  *)>««" 


Figure  5:  Diagram  of  Back  Lighted  Photography 


Figure  8:  Diagram  of  SchJierea  Imaging  System 
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Figure  9:  Unexcited  chamber  pressure,  gas  phase 
velocity  =  93  m/s,  liquid  velocity  »  15  m/s, 
equivalence  ratio  =7.4 
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Figure  iO;  Power  Spectrum,  unexcited,  gas  phase 
locity  =  93  m/s,  liquid  velocity  =  15  m/s 
equivalence  ratio  =7.4  y  '  ' 
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Figure  11:  Excited  chamber  pressure,  gas  phase 
velocity  =  93  m/s,  liquid  velocity  -  15  m/s, 
equivalence  ratio  *  7.4 


Figure  12:  Power  Spectrum,  excited,  gas  phase 
velocity  =  93  m/s,  liquid  velocity  =  15  m/s, 
equivalence  ratio  *  7.4 
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Figure  13:  Unexcited  chamber  pressure,  gas  phas 
velocity  »  93  m/s,  liquid  velocity  *  15  m/s, 
equivalence  ratio  =3.8 


Figure  14:  Power  Spectrum,  unexcited,  gas  phase 
velocity  =  93  m/a,  liquid  velocity  -  15  m/s, 
equivalence  ratio  =3.8 


Figure  15:  Excited  chamber  pressure,  gas  phase 
velocity  -  93  m/s,  liquid  velocity  *  15  m/s, 
equivalence  ratio  -  3.8 


I  5M  1(M  ’$M  ](M  nm  KM  ISM  «M  <SM  MM 

ItWI 


Figure  16:  Power  Spectrum,  excited,  gas  phase 
velocity  »  93  m/s,  liquid  velocity  *  15  m/s, 
equivalence  ratio  -  3.8 


